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Abstract 

This paper is a continuation of a previous one [L. Samaj, J. Stat. 
Phys. 120:125 (2005)] dealing with the renormalization of a guest 
charge immersed in a two-dimensional logarithmic Coulomb gas of 
pointlike ± unit charges, the latter system being in the stability- 
against-collapse regime of reduced inverse temperatures < j3 < 2. In 
the previous work, using a sine-Gordon representation of the Coulomb 
gas, an exact renormalized-charge formula was derived for the special 
case of the pointlike guest charge Q, in its stability regime @\Q\ < 2. 
In the present paper, we extend the renormalized-charge treatment to 
the guest charge with a hard core of radius a, which allows us to go 
beyond the stability border j3\Q\ = 2. In the limit of the hard-core 
radius much smaller than the correlation length of the Coulomb-gas 
species and at a strictly finite temperature, due to the counterion con- 
densation in the extended region (5\Q\ > 2, the renormalized charge 
Q ren turns out to be a periodic function of the bare charge Q with 
period 1. The renormalized charge therefore does not saturate at a 
specific finite value as \Q\ — > oo, but oscillates between two extreme 
values. In the high-temperature Poisson-Boltzmann scaling regime of 
limits (3 — > and Q — > oo with the product (3Q being finite, one re- 
produces correctly the monotonic dependence of (3Q ren on (3Q in the 
guest-charge stability region (3\Q\ < 2 and the Manning-Oosawa type 
of counterion condensation with the uniform saturation of PQ Te n at 
the value 4/7r in the region P\Q\ > 2. 
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1 Introduction 



The concept of renormalized charge is of primary importance in the equi- 
librium statistical mechanics of colloids, see, e.g., refs. [TJ E3 HI El • This 
concept is based on the assumption that, at a finite temperature, the electric 
potential induced by a "guest" (say colloidal) charged particle, immersed in 
an infinite electrolyte, exhibits, at large distances from this particle, basi- 
cally the screening form given by the high-temperature linear Debye-Hiickel 
(DH) theory. From the physical point of view, due to strong electrostatic 
interactions, the guest particle of bare electric charge Q attracts oppositely 
charged electrolyte particles (counterions) to its immediate vicinity, and this 
decorated object may be considered as a new entity of lower renormalized 
charge Q ren . The idea of renormalized charge was confirmed in the frame- 
work of the nonlinear Poisson-Boltzmann (PB) approach there is 
a renormalized-charge prefactor, Q Ten , to the usual Yukawa decay which is 
different from the bare charge Q of the guest particle. An interesting point is 
that, as the absolute value of Q increases to infinity, the renormalized charge 
saturates monotonically at some finite value Q^- The possibility of a more 
general phenomenon of potential saturation was studied from a general point 
of view in ref. [7j and on the exactly solvable 2D Coulomb gas at the Thirring 
free-fermion point in ref. [S|. 

The rigorous validity of the PB approach is restricted, under certain 
conditions making the nonlinear PB theory superior to the linear DH one 
jH], to a specific scaling regime of the infinite-temperature limit [TU]. To 
go beyond this mean-field description, one has to incorporate electrostatic 
correlations among the electrolyte particles, like it was done, e.g., in refs. 
[TTJ H2] • Such approaches always involve some plausible, but not rigorously 
justified, arguments and approximation schemes. 

Another strategy is to concentrate on simplified models which keep the 
Coulomb nature of particle interactions and simultaneously admit the exact 
solution. Specific two-dimensional (2D) Coulomb systems with logarithmic 
pairwise interactions among charged constituents, where the electrolyte is 
modelled by an infinite symmetric Coulomb gas, belong to such category of 
models. The 2D Coulomb gas of ± unit pointlike charges is stable against the 
thermodynamic collapse of positive-negative pairs of charges at high enough 
temperatures, namely for (3 < 2 with (3 being the (dimensionless) inverse 
temperature. In this stability region, the equilibrium statistical mechanics 
of the Coulomb gas is exactly solvable via an equivalence with the integrable 
(l+l)-dimensional sine-Gordon theory; for a short review, see ref. |13j . An 
extension of the exact treatment of the stable 2D Coulomb gas to the pres- 
ence of some pointlike guest charge (s) was done in ref. [Tl]. The problem of 
one (two) guest charge (s) immersed in the 2D Coulomb plasma was shown 
to be related to the evaluation of one-point (two-point) expectation values of 
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the exponential field in the equivalent sine-Gordon model. Based on recent 
progress in the latter topic, two main results were obtained. Firstly, an ex- 
plicit formula for the chemical potential of single guest charge Q was found 
in the guest-charge stability (no collapse of Q with a unit plasma counte- 
rion) region /3\Q\ < 2. Secondly, the asymptotic large-distance behavior of 
an effective interaction between two guest charges Q and Q' was derived. 
As a by-product of this result, considering that Q' corresponds to either +1 
or —1 charged Coulomb-gas species, the concept of renormalized charge was 
confirmed and the explicit dependence of Q Ten on f3 and Q, valid rigorously 
in the whole stability range of \Q\ < 2/(3, was established. For a fixed (3 < 2, 
the renormalized charge Q ren , considered as a function of, say positive, bare 
charge Q, exhibits a maximum at Q — 2/(3— 1/2: this non-monotonic behav- 
ior resembles the one observed in the Monte-Carlo simulations of the salt-free 
(only counterions are present) colloidal cell model JJJ. At the stability bor- 
der Q = 2/(3, Q r en attains a finite value. This fact was an inspiration for a 
conjecture, referred to as "regularization hypothesis", about the possibility 
of an analytic continuation of the formula for Q ren to the collapse region 

Q>2/(3\TQ- 

The validity of the regularization hypothesis was put in doubts by Tellez 
[To] who calculated in detail the formula for the renormalized charge within 
the 2D nonlinear PB theory. The PB approach describes correctly, under 
certain conditions [TJ]j, the scaling regime of limits (3 — > and Q — > oo with 
the product (3Q being finite. As was expected, in the guest-charge stability 
region (3\Q\ < 2, the exact formula for Q ren [JJ], taken with (3 — > 0, was 
reproduced. When (3\Q\ > 2, a hard core of radius o > 0, impenetrable for 
electrolyte particles, has to be attached to the guest charge Q in order to 
prevent its collapse with electrolyte counterions. For the determination of 
the renormalized charge in this regularized case, the connexion problem of 
the PB equation to relate the large-distance behavior of the induced electric 
potential with its short- distance expansion jTB] is of primary importance. 
The numerical PB results of ref. jT3] show that Q ren is always an increasing 
function of the bare Q, and saturates at a finite value as Q — > oo. In par- 
ticular, when the dimensionless positive hard-core radius a = ko {k denotes 
the inverse Debye length) is very small, a — > 0, the renormalized charge 
saturates at the value given by 

(3Q(f e f = - for all values of (3Q > 2. (1.1) 

7T 

This is a manifestation of the Manning- Oosawa counterion condensation ^7J 
ITS] known in the theory of cylindrical polyelectrolytes. 

The study of the renormalization charge within the nonlinear PB ap- 
proach [T3] is trustworthy, however, there are two open problems. Firstly, 
the rigorous validity of the PB theory was proved for 3D electrolytes in 
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the presence of some continuous domains of external charge distributions 
[TUj . The hard-core interaction between the external guest charge and the 
Coulomb-gas particles, which is so relevant in the questionable guest-charge 
collapse region (3\Q\ > 2, was not considered in Kennedy's proof. The sec- 
ond problem, which is probably even more important, is related to the limit 
(3 — > considered in the PB approach. This limit represents a very strong 
restriction which prevents one from seeing nontrivialities in the plot of Q TeD 
versus Q appearing at a strictly finite (nonzero) (3, like the previously men- 
tioned existence of the maximum extreme at the point Q = 2/ (3 — 1/2 [T%] . 
It is evident that the monotonic increase of (3Q rcn as the function of (3Q to 
its saturation value, predicted by the PB theory, is certainly not true for a 
finite inverse temperature (3. 

The aim of the present paper is to extend the exact renormalized- charge 
treatment of ref. to the guest charge with a hard core of radius a, which 
allows us to go beyond the stability border (3\Q\ = 2 of the pointlike guest 
charge. By the spirit, the applied method is similar to the one used in ref. ^H] 
to include hard cores around charged species of the infinite 2D Coulomb gas 
itself. In the limit a — > 0, i.e., when the hard-core radius of the guest particle 
is much smaller than the mean interparticle distance of the Coulomb-gas 
species, due to the counterion condensation in the extended region (3\Q\ > 2, 
the renormalized charge Q ren turns out to be a periodic function of the bare 
charge Q with period 1. The renormalized charge therefore does not saturate 
at a specific finite value as \Q\ — > oo, but oscillates between two extreme 
values. Such behavior indicates that, for a strictly nonzero (3, the Manning- 
Oosawa counterion condensation phenomenon [TT| IT5] should be revisited. In 
the scaling PB regime, one recovers correctly the standard results including 
the saturation formula (jl.lj) . 

The paper is organized as follows. Section 2 reviews the known exact 
information about the equilibrium statistical mechanics of the infinite 2D 
Coulomb gas, including the complete thermodynamics and both short- and 
large-distance asymptotic behaviors of two-point correlation functions. In 
this section, we introduce the notation and present important formulas which 
are extensively used throughout the whole paper. Section 3 deals with the 
chemical potential of the guest charge immersed in the 2D Coulomb gas. The 
case of the pointlike guest charge is analyzed in Section 3.1, the inclusion of 
the hard-core region around the guest charge is the subject of Section 3.2. 
The renormalization of the guest charge is studied in Section 4. In Section 
4.1, the renormalization of the pointlike guest charge is briefly reviewed fol- 
lowing the derivation of ref. [Hj. In Section 4.2, relevant a-corrections due 
to the presence of the hard core are systematically generated in the formula 
for the renormalized charge. A recapitulation is given in Section 5. 



4 



2 Bulk properties of the 2D Coulomb gas 



2.1 Integrability 

We consider a classical (i.e. non quantum) Coulomb gas consisting of two 
species of pointlike particles with opposite unit charges q G {+1,-1}, con- 
strained to an infinite 2D plane A of points r e R 2 . The interaction energy 
of a set of particles {<£,-, r 3 -} is given by J2j<kQjQk v (\ r j ~ T k\), where the 
Coulomb potential v(r) = — ln(|r|/r ) (the free length constant r will be 
set to unity for simplicity) is the regular solution of the 2D Poisson equa- 
tion Av(r) = —2irS(r). The equilibrium statistical mechanics of the system 
is usually treated in the grand canonical ensemble, characterized by the di- 
mensionless inverse temperature (3 and by the couple of particle fugacities 
z + and Z-. Since the length scale r was set to unity, the true dimension of 
z± is [length] -2 +(/ 3 / 2 ). The bulk Coulomb gas is neutral, and thus its ther- 
modynamic properties depend only on the combination ^/z + z_ [20] • It is 
therefore possible to set z + = z_ = z\ however, at some places, in order 
to distinguish between the + and — charges, we shall keep the notation z±. 
The system of pointlike particles is stable against the UV collapse of positive- 
negative pairs of unit charges provided that the corresponding Boltzmann 
weight exp[/3u(r)] = |r| _/3 can be integrated at short distances in 2D, i.e., at 
high enough temperatures such that (3 < 2. In what follows, we shall restrict 
ourselves to this stability range of inverse temperatures. 

The grand partition function 5 of the 2D Coulomb gas can be turned via 
the Hubbard-Stratonovich transformation (see, e.g., ref. [21]) into 

/2tyexp[-S(0)] [ZA) 

with 



S(z) = I d 2 r 

J A 



■(V0) 2 - 2zcos(bcf>) 



l(i;v ' " ' 

being the Euclidean action of the (1 + l)-dimensional sine-Gordon model. 
Here, 0(r) is a real scalar field and / T>(b denotes the functional integra- 
tion over this field. The sine-Gordon coupling constant b depends only on 
the inverse temperature (3 of the Coulomb gas. The fugacity z is renormal- 
ized by the (diverging) self-energy term exp[(3v(0)/2] which disappears from 
statistical relations under the conformal short-distance normalization of the 
exponential fields 

(e^e-""^) ~ |r - r'|" 4b2 as |r - r'| -> 0. (2.3) 

For b 2 < 1 (j3 < 4), the discrete symmetry <fi —>■ <p + 2-nnjb (n being an 
integer) of the sine-Gordon action ()2.2|) is spontaneously broken and therefore 
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the sine-Gordon model is massive [22]- Its particle spectrum consists of one 
soliton-antisoliton pair (S, S) with equal masses M, which coexist in pairs, 
and of S — S bound states, called "breathers" {Bj;j = 1, 2, • • ■ < l/£}, whose 
quantized number at a given b 2 depends on the inverse of the parameter 

The mass of the Bj breather is given by 

m i = 2Msin^yjJ (2.5) 

and this breather disappears from the sine-Gordon particle spectrum just 
when rrij = 2M, i.e., £ = Note that the lightest B\ breather is present 
in the spectrum up to the collapse point 6 2 = 1/2 (/5 = 2). 

The 2D sine-Gordon model is an integrable field theory, so that any mul- 
tiparticle scattering S-matrix factorizes into a product of explicitly available 
two-particle S-matrices satisfying the Yang-Baxter identity [22j. Basic char- 
acteristics of the underlying theory were derived quite recently by using the 
method of Thermodynamic Bethe ansatz. In particular, the (dimensionless) 
specific grand potential u, defined by 

-u = lira -JrlnH(z), (2.6) 

|A|^oo | TV I 

was found in ref. [21]: 

-u= m ) - v (2.7) 
8sin(vr0 1 ' 

Under the conformal normalization of the exponential fields (12.3)) . the rela- 
tionship between the fugacity z and the soliton/antisoliton mass M reads 



r(6 2 



7iT(l - b 2 ) 



A//V /?r((i + o/2) 



2r(e/2) 



2(l-fe 2 ) 

(2.8) 



where T stands for the Gamma function. 



2.2 One-point densities 

The homogeneous number density of the Coulomb-gas species of one sign 
q = ±1, n q (r) = n q with r e R 2 , is defined standardly as the thermal average 
n q = (J2j ^q,qj^( r ~ r j))/3- The charge neutrality of the system implies that 
n + = n_ = n/2, where n denotes the total number density of particles. The 
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species densities are expressible as field averages over the sine-Gordon action 
Q) as follows 



z d(—u) 
2 dz 



±1: 



(2.9) 



This equality and the relation f)2.7j) determine explicitly the density-fugacity 
relationship 



n l-(/3/4) 



2 *1 



P/ A tvi (aiA\\ r /1 a a \i^-W*) 



r(i - (g/4)) 
r(i + (/?/4)) 



I 2' 4-/3' 



2(4-/3) 



;1 



(2.10) 

where F = 2-^1 is the hypergeometric function. Based on this density- 
fugacity relationship, the complete thermodynamics of the 2D Coulomb gas 
was derived in the whole stability regime of pointlike charges (3 < 2 [26J. 

The excess (i.e., over ideal) chemical potential of the Coulomb-gas species 
q = ±1, /ig X , is given by 



exp (-Pfi 



n„ 



), q = ±l. 



(2.11) 



Let fi'g with arbitrarily valued real Q represents an extended definition of 
the excess chemical potential: //q is the reversible work which has to be done 
in order to bring a pointlike guest particle of charge Q from infinity into the 
bulk interior of the considered Coulomb gas. It was shown in ref. |T3] that 
/ig x is expressible in the sine-Gordon format as follows 



exp (-(3^1 



(e iW ) 



(2.12) 



When Q = ±1, one recovers the previous result ()2.11j) valid for the Coulomb- 
gas constituents. Due to the obvious symmetry relation (e ia ^) = (e~ ia< ^) valid 
for any real- valued parameter a, it holds that /ig = H e *n- 

A general formula for the expectation value of the exponential field (e m ^) 
was conjectured by Lukyanov and Zamolodchikov j2H|. In the notation of 
Eq. ()2.12|) . a = Qb, their formula reads 



(e lW ) 



nzT(l - b 2 
f{b 2 



2M {Qb?/(i-b 2 ) 



exp [I h {Q)\ 



(2.13) 



with 



h(Q) 



dt 



sinh 2 (2Q6 2 t) 



2r2 -2t 



t [2 sinh(6 2 t) sinh(t) cosh[(l - b 2 )t] 



2Q z b z e 



(2.14) 
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The integral (l2~T3D is finite provided that \Q\ < l/(2b 2 ): at \Q\ = l/(2b 2 ) the 
integrated function behaves like 1/t for t — > oo what causes the logarithmic 
divergence. In the Coulomb-gas format, the interaction Boltzmann factor of 
the guest Q charge with an opposite unit plasma counterion at distance r, 
r -/3]Q] ) j s integrable at small 2D distances r if and only if /3\Q\ < 2. In terms 
of the sine-Gordon coupling constant b 2 = /3/4, the stability region for \fq 
is therefore expected to be \Q\ < l/(26 2 ) and the couple of Eqs. (|2.13|) and 
(I2.14J1 passes the collapse test. 



2.3 Two-point densities 

At the two-particle statistical level, one introduces the two-body densities 
■ — (J2j^k 3q,qj3( r ~~ r j)^q',qk^( r> ~ r k))i3 which are translationally 



invariant in the infinite 2D space, n^/(r, r') = nS(|r — r'|). The two-body 
density is expressible as an average over the sine-Gordon action (|2.2|) as 
follows 

n$(r) = ^l^ '^); q, q' = ±1. (2.15) 

In close analogy with the previous case of one-body densities, it was shown in 
ref. ^3] that the effective interaction energy of two guest charges immersed 
in the 2D Coulomb gas is expressible in terms of more general two-point 
correlation functions of exponential fields: ^e 1< ^ 6< ^ V < ^' b ^ (r ' ) ) with real- valued 
charge parameters Q and Q' . A systematic generation of the short- and large- 
distance asymptotic expansions for these two-point correlation functions is 
available with the aid of special field-theoretical methods. 

The short-distance expansion of (e 1 ^ 6 ^ -^ 1 ^ b( ^ r )) can be obtained by using 
the method of Operator Product Expansion (OPE) [27 . The OPE has the 
form 



e i W (o) eW (r) = £ {c%°,(r)eW +Q ' +n W(0) + •••}, (2.16) 

n=— oo 

where the dots stand for subleading contributions of the descendants of 
e i(Q+Q'+n)^ like (00)2^2 e i(Q+<y+n)ty eta The coefficients C read 

Cg'g,(r) = /Z l n lr 4f>2[QQ ' +n(Q+Q ' )+n2/2]+2|n|(1_fe2) /3Q/ (VV 4 " 4fe2 ) , (2.17) 
where the functions / admit analytic series expansions 



oo 

/aQ'(*)=£/; ,0 (3>w- (2.i8) 

Each coefficient /"'° is expressible as a 2(\n\ + j)-fold Coulomb integral de- 
fined in the infinite plane R 2 . The leading terms fo'°{Q,Q') in the series 
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(|2.18|) are expressible as 

C(Q,Q') = 1, 

fo'°(Q,Q') = UQb 2 ,Q'b 2 ,b 2 ) forn>0, (2.19) 
»[( 

Here, 



fo'°(Q,Q') = J\n\(-Qb 2 ,-Q'b 2 ,b 2 ) forn<0. 



1 /■ n 

Jn (a,a',b 2 ) = -JU (d 2 r>/ a |l - r/°') [] |r, - r fc | 4fe2 (2.20) 



j=i j<k 



with 1 being a point on the unit circle, say (1, 0). This integral is convergent 
if and only if the parameters (a,a',b 2 ) fulfill the inequalities a > —1/2, 
a' > -1/2 and a + a' < -(n - l)b 2 - 1/2. The integral (j2.20|) was evaluated 
by Dotsenko and Fateev j3UJ |3T] : 



j n (a } a',b ) 



n-l 

' 1,2 



71 



II 7(J& 2 ) II 7(1 + 2a + A;6 2 ) 7 (l + 2a' + kb 2 

j=l k=0 



lip 2 ) 

x 7 (-1 - 2a - 2a - (n - 1 + k)b 2 ) . (2.21) 



Hereinafter, we use the notation j(t) = T(t)/T(l — t). The result (|2.21j) 
represents an analytic continuation of the integral ()2.20|) to all values of the 
parameters (a,a',b 2 ). We would like to emphasize that the OPE algebra 
(I2.16|) is the operation which can be used in any multi-point correlation 
function of exponential fields to reduce its order as soon as a couple of points 
is close to one another. 

The large-distance asymptotic expansion of the truncated two-point cor- 
relation functions 

^ e iQ&<«0) e iQ'6<«r)^T = ^iQb^^Q'b^r)^ _ ,£QM>\ /gQ'ty) (2.22) 

can be obtained by using the form-factor method (221- The form-factor rep- 
resentation is formally expressed as an infinite convergent series over multi- 
particle intermediate states, 

00 1 x . r°° d0i . . . d0/ 



iV! eit ... ;eN J-00 {ITT) 

l F Q ,(6 N , . . . , 9 1 )e~ r ^ m ^ cosh9 >. (2.23) 



N=l ' ei,---,ejv ' 
ejV,...,e a 



Here, e indexes the particles of the sine-Gordon spectrum (say e = +/— for a 
soliton/antisoliton and e = j for the Bj breather), the rapidity 6 G (—00, 00) 
parameterizes the energy and the momentum of the particles, and F denotes 
the corresponding form factors. In the large-distance limit r — * 00, the 
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dominant term on the rhs of Eq. (|2.23|) corresponds to the intermediate state 
with the minimum value of the total particle mass X/JLi m ej , at the point 
of vanishing rapidities. In the stability region of the pointlike Coulomb gas 
< b 2 < 1/2, the lightest particle, which can exist in the spectrum alone, is 
the B\ breather of mass m\. For this particle, the one-particle form factors 
F Q (9)i and 1 F Q ,{6) = F QI (6) 1 were calculated in refs. [33] and [H] : 

F Q {6) 1 = -i(^)y/^^ 1 (2.24) 

sm«) 



A 



where 

- sin«) cos (Hf\ exp f- f* **t\ 
7T \ 2 / \ Jo 7T sm t J 

and £ is defined in Eq. (|2.4|) . Since the form factor (|2.24|) does not depend 
on the rapidity, the integration over 6 in ()2.23|) can be done explicitly by 
using the relation 

r°o c ]f) / it \ V 2 

Lt 6 "™' ~" = *-<•»"■> —G^;) e ' m "- (2 - 26 > 

Here, K is the modified Bessel function of second kind. Consequently, 

/ AQb<fi(0) iQ'6<Kr)\T / ~ \ 1/2 

K , ' [Q}[Q'}\{ — — ) e~ mir , (2.27) 

where the symbol [Q] stands for the ratio 

[0] = ^g>. (2.28) 
sm(vr£) 

We see that, at large distance r, the truncated two-point correlation func- 
tion factorizes into the product of separate phase contributions [Q] and [Q'\. 
The inverse correlation length in the exponential decay, mi, is determined 
exclusively by the Coulomb-gas system. Using Eqs. (j2.5|) - ()2.10|) . mi is 
expressible as 

[ gin (jm zM ] ' /2 (2 29) 

Here, k = ^/2-Kf3n denotes the inverse Debye length; in the high-temperature 
(5 — > limit one has mx ~ k. The formula ()2.29|) describes the renormaliza- 
tion of the inverse correlation length at a finite temperature. 



10 



3 Chemical potential of the guest charge 



3.1 Pointlike guest charge 



The chemical potential of the pointlike, without any loss of generality say 
positive, guest charge Q > 0, immersed in the bulk of the 2D Coulomb gas, 
is given by the relation ([2.12)1 . complemented by the explicit formulas ()2.13|) 
and (|2.14|) for the expectation values of exponential fields in the sine-Gordon 
model. Eqs. ()2.13j) and (|2.14|) are valid provided that Q < Q c , where 



Qc 2l)2 




(3.1) 



is the "collapse" value of the guest charge. Note that Q c > 1 in the considered 
stability regime of the Coulomb gas (3 < 2 (6 2 < 1/2). 

Eqs. (I2.13j) and (|2.14j) can be continued analytically to the region Q > Q c 
by using the reflection method developed in refs. j^HlEnj- Exploring a close 
relationship between the Liouville and sine-Gordon theories, it was argued 
that the sine-Gordon expectation values (e 10 ^) obey the relation 



<e ia ^> = J R(a)(e- i(a+w ), 
where Q = 6 _1 — b, and the reflection amplitude reads 



(3.2) 



R(a) 



7izT(l - b 2 ) 
b 2 T(l + b 2 ) 



(2a+g)/b rfl + ^ + |)r(l- 2ab - Qb) 



2ab + Qb) 



r i 



!-! r( 



(3.3) 



Eqs. ([2.13)1 and ()2.14j) . in the notation of a = Qb restricted to \a\ < 1/(26), 
represent just the "minimal" solution of the reflection relations ()3.2)1 and 
()3.3j) . Using successively the reflection formula ()3.2j) . any expectation value 
(e 101 ^) with \a\ > 1/(26) can be reduced to a product of amplitudes of type 
()3.3j) and the expectation value of the exponential field with a phase whose 
absolute value < 1/(26). 

For our purpose, we take advantage of the symmetry property (e m< ^) = 
(e~ m< ^) and rewrite the reflection relation ()3.2j) . taken with a = Qb, in the 
following way 



) = ( 



-iQb<t>\ 



e 

i?(-Q6)(e i(Qfe - g) *>, 



where 
R(-Qb) 



7CZT{1 - b 2 ) 

6 2 r(i + b 2 ) 



-i 2Q+1-6- 



r(-2Q + 6- 2 )r(2Q6 2 + 6 2 ) 
T(2 + 2Q- 6" 2 )r(2 - 2Q6 2 - 6 2 ) 



(3.4) 



(3.5) 
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While the minimal relations (|2.13j) and (|2.14|) are applicable to — Q c < Q < 
Q c , the "first-order" relations (JO) and dSZHJ) hold for -Q c <Q-{Q/b)< Q c , 
which is equivalent to 

Q c - 1< Q < Q c + (b- 2 - 1). (3.6) 

In the considered stability region b 2 < 1/2, one has Q c > 1, so that the 
lower bound for the guest charge Q in Eq. (|3.6j) is a positive number, and 
b~ 2 — 1 > 1, so that the upper bound for Q in Eq. (|3.6j) is larger than Q c + 1. 

The only singularities (in particular, simple poles) of the reflection am- 
plitude R(—Qb) in Eq. ()3.5|) come from the Gamma function r(— 2Q + b~ 2 ), 
when its argument attains zero or a negative integer, i.e., at the infinite 
sequence of equidistant Q-points 

ii 

Q* n = Q c + -; n = 0,l,2,.... (3.7) 

The first three singular Q-points, which will be of special interest in what 
follows, explicitly read 

Q*o = Qc: Ql = Qc + \, Q* 2 = Q c + l; (3.8) 

note that these points lie in the range ()3.6|) of the applicability of the first- 
order reflection relations ()3.4|) and (|3.5j) . Using the residuum formula for the 
Gamma function 

(-l) n 1 

lkr-n + e=^- n = 0,1,2,...), 3.9 

Eqs. ()3.4|) and (|3.5j) imply, as e — >• + , the following singular behaviors 

(JQH\ „ H_ /j(Q* -i)b<t> 



) ~ w- e \ e / ioi Q = Qo-^ (3-10) 



/ 4e 



rf(l - b 2 ) 
b 2 T{l + b 2 ) 



T(l + 2b 2 ) 
T(l -26 2 ) 



x ^(Ql-^ for g = Q* _ C) (3.11) 

and so on. 



3.2 Guest charge with hard core 

Let us introduce an impenetrable hard disc of radius a around the positively 
charged Q > guest particle localized in the bulk interior of the 2D Coulomb 
gas, say at the origin r = 0. The radius a has the dimension of length; the 
most natural dimensionless quantity involving a is chosen as 

a = rriicr, (3-12) 
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where m x is the inverse correlation length of the Coulomb-gas constituents 
given by Eq. ([2.290 . The corresponding excess chemical potential of the 
guest charge will be denoted by yUg (a). As soon as o > 0, yUg ( a ) m ust be 
finite for any value of Q. It was shown in ref. ^1] that ^ < q(<j) is expressible 
in the sine- Gordon format as follows 



where the average (• 

<•">* 
S a {z) 



exp \ -(3jJinirf \ 

is defined by 
1 



H(z) 

+ z 



V(j)e~ Sa{z) 
d 2 r 



r<cr 



(r) 



(3.13) 

(3.14) 
(3.15) 



Here, S(z) is the usual sine-Gordon action (j2.2jl . 

The averaged quantity (e l( ^ b ^°^) a can be formally expanded around S^z): 



n=l 



(3.16) 



where 




J$>M = 


/ 




Jr<a 


4 2 V) = 


J 

Jr<a 







d 2 r 
d 2 r 



/ e iy60(U; e i60(.rj\ _|_ ^ e iy60(Uj e -i6^)(rj\ 



dV 



r'<a 
e -i6< 



/ e iQb?i(O) e ife0(r) e ife0(r')^ _j_ / e iy6<A(,Oj e i6^(rj e -i6</>(r'j\ 



(r) o i60(r')\ 



3 iQ6<K0), 



(3.17) 

iQ^O)^,/.. „:,•), ,-i5<«r')\ 
i6^(r')\" 



(3.18) 



and so on. The expansion (J3.16[) represents the basis for a systematic gen- 
eration of a-corrections to the pure sine-Gordon expectation value (e 1 ^^), 
in the limit a — > 0. To be more particular, for a given value of the guest 
charge Q, the rhs of Eq. (|3.16[) contains a-dependent terms which either 
vanish (and so they are irrelevant and omitted) or diverge (and so they are 
relevant and preserved) in the considered limit a —>■ 0. Since the chemi- 
cal potential /ig(cr > 0) in Eq. (j3.13j) is finite for any Q, these relevant 
a-dependent terms must cancel the singularities of (e 1 ^ 6 ^) occurring at the 
points {Q*}^1 given by Eq. (|3.7|) . The main advantage of the expansion 
(I3.16[) is that, as Q increases from to oo, the singularities of (e 1( ^) at 
the points Qq, Q\, . . . are eliminated successively via the respective integral 
., where the multi-point correlation functions are evaluated 



(1) r (2) 



terms I { 

by using the short-distance OPE scheme described by Eqs. (j2.16[) - ([2.21l) . 
We shall document this claim on lower levels in the next paragraphs. 

In the stability region of the pointlike guest charge < Q < Qq, one has 



,iW(0) 



} (j ~ <e iW > 



as a 







(3.19) 
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in the sense that all a-dependent terms on the rhs of Eq. (|3.16|) vanish in 
the limit a — > 0. This is a direct consequence of the fact that the chemical 
potential /ig(cr), given by Eq. ()3.13|h is finite even for a = in the guest- 
charge stability region \Q\ < Qq. 

When Qq — e < Q < Ql (in what follows, under e we shall understand 
a positive number going to zero, e — ► + ), the first term on the rhs of Eq. 
(|3.16|) has to be taken into account. The corresponding integral Iq (er), 
defined in Eq. (I3.17|) . consists of two parts: the first term corresponds to 
the correlation function of the pointlike guest charge Q with one coion (unit 
Coulomb- gas charge of the same sign as Q), the second term corresponds 
to the correlation function of the Q-charge with one counterion. In the 
considered limit a — > 0, the systematic short- distance expansions of the 
underlying correlation functions, based on the OPE scheme (|2.16)) - (j2.21j) . 
has to be applied. The dominant contribution comes from the correlation 
function of the Q charge and one counterion, in the leading short- distance 
order corresponding to the bare-Coulomb Boltzmann factor: 



e iW(0) e -i6<«rh ^ / e i(Q-l)60\ r -4Q6^ (3.20) 



Consequently, 



, i(Q-l)6<M 2.2-4Q6 2 

I«(a)~ 27 r ^_ 4 ^ 2 _ 4g&2 asa-0. (3.21) 



After simple algebraic manipulations, Eq. (|3.16|) then reads 



/ i(Q-l)6<M fr4b 2 (Q*-Q) 



I iQtxt>(0)\ _ _S1 LI (3 92) 

{e )(7 ^o {C 1 26 2 m fW5-Q) Ql-Q ■ [ ] 

Close to the collapse value of the guest charge, i.e., when Q = Qq — e, one 
expands 

-4b 2 (Q* -Q) 1 

— — - = -+46 2 lna + 0(e). (3.23) 

VO — V e 

Inserting this expansion into Eq. (|3.22j) . the leading singular term of order 
e _1 cancels exactly with its counterpart in (e lQ6< ^) [see Eq. (|3.10|) ]. and one 
ends up with the logarithmic dependence on the short-distance cutoff a: 

(e i( 3S^(o)) a 2irz(e i{Q o- 1)brp )\na. (3.24) 

The outlined cancellation of singularities at Qq, due to the presence of the 
hard core around the guest charge, is an important check of the consistency 
of the present expansion procedure. For Q > Qq, the exponent of a in Eq. 
f!3.22|) is negative as was anticipated. 
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When Ql — e < Q < Q* 2 , in the limit a — ► 0, two additional relevant 
er-contributions arise on the rhs of Eq. (|3.16|) . The first contribution has the 
origin, as in the above paragraph, in the integral (|3.17|) as the result 

of the next-to-leading expansion term for the correlation function of the Q 
charge with one counterion: 



/ e iW(0) e -i6*(r)v _ /eW-D^^r" 4 ^ 2 + ^' 1 (-g6 2 ,6 2 ,6 2 )(e i ( Q - 2 ^)r 2+4fe2 - 8Q62 . 

r— >0 

(3.25) 

This expansion was derived from the OPE scheme ()2.16|) - (J2.21J) assuming 
that Q > 1, which is indeed true for the considered range of Q values. Thus, 

/ i(Q-l)6<M -2-4Q6 2 



/ e i(Q-2)b<t>\ frA+Ab 2 -&Qb 2 

X m 4 + 4 b2 -s^ 4 + 462 _ 8 g 62 as --°- ( 3 - 26 ) 

According to the Dotsenko-Fateev result 1)2.21)1 . it holds 

U-Qb 2 , b 2 , b 2 ) = vr 7 (l - 2Q6 2 ) 7 (1 + 26 2 ) 7 (-l - 2b 2 + 2Qb 2 ). (3.27) 

The second relevant a-contribution has the origin in the integral Iq (a) 
f)3~TKj) . namely, in the correlation function ( e iW(o) e -i60(r) e -i60(r')) of the Q 
charge with two electrolyte counterions. At short distances r and r', this 
three-point correlation function is governed by the bare-Coulomb Boltzmann 
factor and behaves like 

/ e iQ6^(0) e -i6^(r) e -i6^(r')\ ^ / e KQ-2)6<M ^Qfe 2 / A -4Q6 2 | r _ r /|4ft 2 _ /g 2g x 

Thus, 

ffVWff { flZlW +4b2 ~ 8Qb2 (3.29) 

(2) 

where Jq is the Coulomb integral defined inside the unit disc: 

4 2) =/ d 2 r / dVr- 4Qfe2 (r')- 4Qfe2 |r-r'| 4fc2 . (3.30) 

Jr<l Jr'<l 

The derivation of its series representation 

-^2 oo i 

.(2) = 7T ^ 1 

Q 1 + 6 2 - 2Q6 2 ^ 1 - 2Q6 2 + j 



r(j - 26 2 ) 
j!r(-26 2 ) 



(3.31) 



is outlined in the Appendix, see the final Eq. (jA.7j) . Inserting the formula 
(13.26)1 together with Eqs. (j3.29j) and ()3.31j) into the basic series expansion 
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(|3.16j) . after simple algebra one finally arrives at 

7TZ (e^O" 1 )^) cr^ 2 (Q*-Q) 



jQb<f>(0) 



)a 



(e lW ) - 



2V m f(Qo-Q) Ql-Q 
i2){ nz) 2 (e 1 ^- 2 )^) a^(Qi-Q) 
Ab 2 m 8b2 (Q*i-Q) Q|-Q ' 



where 



f (2) _ 



i=o 



-2Q6 2 +j 



r(j - 2& 2 ) 
j!r(-26 2 ) 



- -ji(-QV 

7T 



b\b 2 ). 



(3.32) 



(3.33) 



In addition to the previous result ()3.22|) . Eq. ()3.32|) contains the a-dependent 
term whose exponent is negative, and therefore relevant, just for Q > Q\. 
This term removes the singularity of the pointlike expectation value (e 1 ^ 6 ^) 
at Q = Ql described by the relation (j3.11|) . To show this fact, we first recall 
the formula (jA.ll)) derived in the Appendix: 

2 



j=0 

Thus, at Q 



1 



rq - 2b 2 ) 

j\T{-2b 2 ) 



1 

2^ 



(3.34) 



Q*, Eq. ((QH)) yields 



(2) 



2tt 



1 

fe 2 



Q!& 2 ,& 
2 



7 (-6 2 ) 7 (1 + 26 2 ) 



r(i-^) 
r(i + 6 2 ) 



r(i + 2& 2 ) 
r(i -26 2 )' 



(3.35) 



Here, we have used the Dotsenko-Fateev result ()3.27|) and the standard re- 
lation T(a; + 1) = xT(x) for the Gamma functions. When Q = Ql — e, after 
inserting the expansion 



(3.36) 



__ = _ + 84WO(e) 

and the explicit form of /q* ()3.35j) into Eq. ()3.32j) . the leading singular term 
of order e _1 cancels exactly with its counterpart in (e 1 *^) [see Eq. fl3.11j) ]. 
As before at the collapse value Qq [see Eq. f!3.24|) ]. one ends up with a 
logarithmic dependence on the short-distance cutoff a at the point Q*. 

Note that the exponents of the two a-dependent terms on the rhs of Eq. 
(EP2J) satisfy the inequality 46 2 (Q*-Q) < 8b 2 (Ql-Q) up to Q = 2Q|-Q* = 
Q\. This means that, in the a — > limit, the first a-dependent term with 
the exponent 46 2 (Qq — Q) is dominant in the whole interval Qq < Q < Q^: 



jQb4>(0) 



HZ 



2b 2 m f(Qo~Q) Q-Ql 



Qo<Q< Ql 



(3.37) 
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Crossing the next singular point this term becomes subleading and the 
second a-dependent term with the exponent &b 2 (Q\ — Q) takes the dominant 
role, it turns out that up to Q\. 

One can proceed along the above lines to higher-order ^-contributions, 
with an increasing amount of algebraic laboriousness. In this paragraph, we 
indicate the general structure of the relevant a-contributions on the rhs of Eq. 
(J3.16j) . From among correlation functions in the integral Iq (<t), the impor- 
tant one corresponds to the configuration of one guest charge Q surrounded 
by n counterions, (e lQb< K°)e~ l6( M ri ) _ _ _ e -ib<t>(r n )y ^ e leading short- distance 
expansion of this correlation function is governed by the bare-Coulomb Boltz- 
mann factor: 

/ e iQ60(O) e -i60(n) . . . e -ib</>(r„)\ ^ uAQ-n)bM r -iQb 2 _ _ _ r ~4Qb 2 TT | r _ _ r 1 4b 2 

n,...,r n ->0 " *■}-. % 3 

(3.38) 

Rescaling the 2D spatial coordinates r 1; . . . , r n by a and regarding the equal- 
ity 2n - AnQb 2 + 2n(n - l)b 2 = 4n5 2 (Q*_ 1 - Q), one obtains 



^V)~4 n) tL n *\ ^ {Q ^~ Q) as*-0, (3.39) 

where Jq^ is a 2n-fold Coulomb integral defined in the unit disc. This a- 
contribution has its counterparts, which possess the structure of Eq. ()3.39|) 
except of prefactors different from Jq , in each of the lower-order integrals 

Iq\<j), . . . ,Iq 1 ' ) (cr) in the series f!3.16|) due to the existence of higher-order 
terms of the short- distance expansions of the correlation functions under inte- 
gration. In analogy with Eq. (J3.32j) . the summation over all cx-contributions 
belonging to the family (J3.39J) can be formally represented as 

(n) (-2^r (e^-*)*) g^Un-^-Q) 
Jq n!(4n6 2 ) m ^ (Q: _ 1 -Q) Q* n l - Q ' [6 ' W) 

The term fl3.4()|) becomes relevant for Q > Q^-i, the prefactor is such 
that it cancels exactly the singularity of the pointlike mean value (e 1 ^ 6 ^) 
at Qn-i- Two consecutive nth and (n + l)st cr-contributions of type ()3.4U|) 
interchange their dominant role when the respective exponents Anb 2 {Q* n _ 1 — 
Q) and A(n+l)b 2 (Q* n -Q) coincide, i.e., at the point Q = (n+l)Q* -nQl_ ± = 
Q\ n - This means that, in the a — > limit, (e 1 ®^ ^) a is dominated by the 
nth a-contribution (|3.4U|) in the range Q^in-x) — Q < Q\n- 
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4 Renormalized charge 
4.1 Pointlike guest charge 

The pointlike guest particle of charge Q > 0, localized at the origin and 
surrounded by the infinite 2D Coulomb gas, evokes position-dependent den- 
sity profiles n g (r) of the electrolyte species q = ±1. In the sine-Gordon field 
representation, these density profiles are given by [Tl] 

( e iQb<t>(0) iqb0{r)\ 

At large distances from the guest charge, the form-factor result ()2.27|) implies 

/ 7T \V2 

n q (r) - n q ~ -n ? [g][(2]A ^ — - j exp(-mir) as r — > oo. (4.2) 

Here, [g] = g for g = ±1, the function [Q] (|2.28jl is expressible in terms of 
the inverse temperature (3 as follows 

_ S jn[7r J gg/(4 - gj 

M " smM/(4-^)] ' l4 ' dj 
and the parameter A ()2.25|) as follows 

Using Eq. (|4.2jl . the induced charge density p of electrolyte particles, 
defined by p(r) = n+(r) — n_(r), reads 

/ 7T \ V2 

p(r) ~ -n[Q]X exp(-mir). (4.5) 

r->oo \2mirJ 

The average electrostatic potential ■?/> is related to the charge-density profile 
through the 2D Poisson equation, 

A^(r) = -2vrp(r). (4.6) 

Considering the circularly symmetric Laplacian A = r _1 d r (rd r ), the asymp- 
totic formula ()4.5|) gives 

M{r)^jQ]\ (^) 1/2 exp(-m 1 r). (4.7) 



■r— >oo 



In the Debye-Hiickel limit (5 — > 0, it holds [Q] ~ Q, J7ii ~ ft and A ~ (3. Eq. 
(|4.7jl then reduces to 

AMO ~ /3Qf^) V2 exp(- K r). (4.8) 



Eq. (|4.7|) is consistent with Eq. (|4.8|) . in terms of the dimensionless com- 
binations m x r and nr, respectively, when one introduces the renormalized 
charge Q Ten as follows [Hj 

PQr C n=[Q}\( — ) 2 ■ (4.9) 

\miJ 

We see that the dependence of the renormalized charge Q ren on the bare 
charge Q enters only via the function [Q] defined by Eq. ([4.3)1 . The rigorous 
validity of the result ()4.9|) is restricted to the guest-charge stability region 
< Q < Q c - Within this range of Q values, the function [Q] exhibits the 
maximum at Q = Q c — 1/2. It can be readily checked that at the collapse 
value of the guest charge this function satisfies the important relation 

[Q c ] = [Qc-i] (4.10) 

which was not noted in the previous work [Tl]. The relation f)4. lOj) reflects 
the effect of the counterion condensation: the renormalized charge associated 
with the collapse value Q c of the bare charge is identical to that of the 
bare charge Q c — 1 because of the condensation of one counterion from the 
electrolyte onto the guest particle. 



4.2 Guest charge with hard core 

We now consider the positive guest charge Q possessing the hard core of 
radius a. The evoked density profiles n g (r) of the electrolyte species q = ±1 
vanish inside the hard disc, n q (r) = for r < a. Using the formalism of ref. 
[T4"j it can be shown that outside the hard-core region the number densities 
of the electrolyte species are given by the following generalization of the 
relation (ETTD 

/ e iQ6</»(O) e igfe0(r)\ 

n ^) = z forr>a, (4.11) 

where the average (• • -) a is defined in Eqs. (|3.14|) and (|3.15jl . The relation 
(|4.11|) can be rewritten in a more convenient form: 

n q {v)-n q = n q ( e iQ^o)) g ( e i^) ' r > °' (412) 

The one-point quantity (e 1 ^ 9 ^ )^ was analyzed in Section 3.2, the large- 
distance analysis of the two-point correlation functions (e'^^e 196 ^ 1 ")) a (q = 
±1) will proceed analogously. Expanding formally (■ ■ -) a around the sine- 
Gordon action S(z), one obtains 

/ e iQ We i#(r)j ff = (eW (o) e i^(r) ) + y L_!L/W( r; a ) ; (4.13) 
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where 



lW(r;a) = f d 2 n £ ( e ^^(o) e i ?1 ^( ri ) e W(r)^ (414) 

J ri«r qi= ±i 

r; a) = / d 2 n / d 2 r 2 £ E 



r(2) 
'Q,9 



<?i=±l 92= 

/ e iQ6^(0) e i ? i6^(r 1 ) e i ?2 6^(r 2 ) e i96^(r)\ (4-15) 



and so on. As before, increasing the value of the guest charge Q from to 
oo, the expansion (j4.13|) represents the basis for the systematic generation 
of (j-corrections which are relevant (i.e., have a negative exponent) in the 
(7 — ► limit. 

In the stability region of the pointlike guest charge < Q < Qq, one has 



e 



Inserting this relation together with the previously derived one (j3.19|) into 
Eq. (j4.12j) and repeating the procedure of Section 4.1, one ends up with the 
renormalized-charge formula (|4.9jl . 

When Qq— e < Q < Q\, the relevant contribution to the integral Iq\{t; a) 
(I4.14|) comes from the q% = — 1 counterion correlation function. Since the 
counterion position vector r x is close to the origin 0, we can apply the short- 
distance OPE formula 



,iQty(O) e -ib0(n) _ r -4Qfe 2 e i(Q-l)60(O) ) ^jjX 

n— »o 



which implies 



^&(^)~ ^ m 2-^ "^g^ as ^°- ( 4 - 18 ) 

The consideration of this formula in the basic expansion 1)4.13)1 leads to 

7V7 / i(Q-l)&0(O) p igb(£(r)\ ±.46 2 (Q*-Q) 
/ iQty(O) i 9 60(rh / iQ6^(0) W(r)v _ e / ° 

(4.19) 

Close to the collapse point, Q = Qq — e, using the large-distance form- 
factor representation of two-point correlation functions (J2.27)) and the for- 
mula (|4.1Uj) . singular terms of order e _1 disappear from the rhs of Eq. (|4.19|) 
in the same way as it was in the case of (e 1 ^ 6 ^ ))^ see Eqs. (|3.22)) - (|3.24)) . 
Finally, inserting ()4.19|) into Eq. (|4.12j) and using the form-factor represen- 
tation (l2~27l) together with the result (l3~2~2l for (e i<w(0) ) ff , one obtains 

1 /2 

n q (r) -n qr ^-n q [q]g^\a)\(^—^ exp(-m 1 r). (4.20) 
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Here, 



with 



(4.21) 



<W ~ ^)lQ]-w m f^-l cr-Q W- 1 '- (4 ' 22) 

Dg) W ~ -c-*) " < e ' W ' W > ^ W " g) (4 23) 

The denominator Dq 1 \(j) is equal to the one-point average (e 1 ^ 6 ^ )) „ given 
by Eq. ()3.22|) and the numerator Ng\a) possesses the same structure as 

Dq ' (a) but each term is multiplied by the function [■ • ■] ()2.28|) with the 
argument related to the phase of the corresponding averaged exponential 
field; we shall see that this rule takes place also on higher levels. At the 
collapse point Q = Qq, the equality ()4.1()j) implies that ggl (a) ~6-_ > o[(3o], 
and one recovers the pointlike result (|4.2j) . This means that Q ren is the 
continuous function of the bare charge Q at the collapse point Qq. For 
Qq < Q < Ql, the term of order ( j 462 ( < 2o -< 2) is dominant in the limit a — > 
in both Dq (a) and Nq\<t), so that Qq = [Q — 1]. The repetition of the 
procedure of Section 4.1, starting from Eq. ()4.20|) . then leads to 

/3Qrcn = [Q-l]\(— )\ Ql<Q<Q\. (4.24) 
\miJ 

Comparing this relation with Eq. (j4.9|) (valid for Q < Qq) we see that, in 
the considered range of Q values, the bare charge Q of the guest particle is 
effectively reduced by 1 due to the condensation of one electrolyte counterion. 

In the region Ql < Q < Q^, one recovers the formula of type ()4.2())1 with 
the substitution 5^(0") — > g Q 2 \a), where Qq\<j) is again the ratio of type 



flOIJ): 



,(2)^ _ N Q 



" ] = » (425) 



Here, the denominator Dq 2 \<j) is equal to the one-point average (e 1 ^ 6 ^ ^ 
given by Eq. (|3.32|) and, as before, the numerator 



jr-z / p i(Q-l)6</>\ i-4^(QS-Q) 



+fQ Ah 2 m f(Qi-Q) Ql-Q [Q ~ 2] (426) 
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arises from Dq (a) by multiplying each term by the function [■ • •] with the 
argument related to the phase of the corresponding averaged exponential 
field. In the a — > limit, once again the term of order <j 4b2 (Qo _( 2) is dominant 

(2) C2l 

111 both D K Q>{a) and N$' (a) in the considered range Q\ < Q < Q* 2 , so that 

fty\ —————— 

9q = [Q — 1]. The previously derived formula (|4.24j) thus applies to the 
whole region Q* Q < Q < Q 2 , 

/3Qrcn = [Q-1]X(— )\ Q* <Q<Q* 2 . (4.27) 
\mi/ 

It can be shown that the next term in Eq. ([4.26)1 . of order cf 8b2 (Q*-Q) , 
becomes dominant in the region Q 2 < Q < Q\. Consequently, in the limit 
a -> 0, 

/3Q ren = [Q-2]\(— )\ Q* 2 <Q<Qt (4.28) 

This relation describes the simultaneous condensation of two electrolyte 
counterions onto the guest charge Q. 

With the aid of the general analysis outlined at the end of Section 3.2, 
one can extend the treatment, performed above in the range Qq < Q < Q 2 , 
to an arbitrary interval Q^n-i) — Q < Q%.n ( n = 1>2. . .). For a given n, 
the dominant term (in the a — » limit) in (e 1 ^ 1 ^ ))^ is of type ()3.40|) . The 
corresponding factor generated by ^e 1 ^"™-""^) is [Q — n), which implies 



PQren=[Q-n]\^ 2 , Q^n-l) <Q < Q*2n (4-29) 



for n = 1,2,.... The bare charge Q of the guest particle is therefore ef- 
fectively reduced by n due to the condensation of n electrolyte counterions. 
Since Q 2n = Q c + n, the formula ()4.29|) tells us that beyond the collapse 
point Qc, i.e., when f3Q > 2, the renormalized charge Q ren is the periodic 
function of the bare charge Q with period 1. The minimum value of Qren, 
given by 

WLT> = A Bin .Iff-/]^ (4.30) 

is acquired at the infinite sequence of points Qq,Q 2 ,QI, ■ ■ ■', the maximum 
value of Qren, given by 

sm [ir (3 /(A- (3)] 

is acquired at the infinite sequence of points Ql, Q%, Ql, .... It is obvious that, 
in the considered limit a — > 0, the renormalized charge does not saturate at 
a specific finite value as Q — > oo, but oscillates between the two extremes 
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given by Eqs. (ET3UJ) and (IOT|l . Note that Q^ n) is positive, i.e., has the 
sign of the bare Q, in the considered Coulomb-gas stability region < /3 < 2; 
consequently, there is no charge inversion in the system. 

Like for example, the plot of Q rcn versus Q at the inverse temperature j3 = 
1 is presented in Fig. 1. The solid line represents the true dependence given 
by Eqs. (|4.9|) and (|4.29|) . The dashed line corresponds to the "naive" analytic 
continuation of the pointlike formula (J4.9j) into the counterion-condensation 
region Q > 2. 



1 




Figure 1: Qrcn versus the bare charge Q at j3 — 1 (solid line). 

We now check the obtained exact results on the high-temperature PB 
scaling regime of limits (3 — > and Q — > oo with the product (5Q being 
finite. In the guest-charge stability region < (3Q < 2, the consideration of 
mi ~ k, X ~ (3 and [Q] ~ 4[sin(7r/5Q/4)]/(7r/3) in Eq. (gU) results in 

(3Qren = ^ sin {^^j , < /3Q < 2. (4.32) 

Evidently, f3Q reri is an increasing function of (3Q in the stability interval. It 
acquires its maximum value, equal to 4/tt, when f3Q approaches the collapse 
value 2. For /3Q > 2, the two extreme values (j4.3(jp and (j4.31|) coincide in 
the PB scaling regime: 

PQl^ = f3Q^ = -, f]Q>2. (4.33) 
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This is equivalent to saying that Q ren saturates at the value given by PQ^l^ = 
4/ V for all (3Q > 2, which is in agreement with the result (jl.lj) of the 
Manning-Oosawa theory of counterion condensation [15J. We would like to 
emphasize that both the monotonic increase of /3Q ren in the stability region 
< (3Q < 2 and the saturation of f3Q ren at the uniform value for j3Q > 2 
are the specific features of the f3 — > limit. In the case of a strictly fi- 
nite (nonzero) value of the inverse temperature /3, the previously described 
non-monotonic oscillatory dependence of Q ien on the bare charge Q takes 
place. 

5 Conclusion 

Exact results of ref. ^3], which concern thermodynamic characteristics of 
a pointlike guest charge Q (say, > 0) immersed in the infinite 2D Coulomb 
gas of ±1 charged species, are rigorously valid in the stability region of the 
guest particle < (3Q < 2. The analytic continuation of these results to 
the "collapse" region (3Q > 2 turns out to be inadequate. In the region 
(3Q > 2, in order to prevent the direct collapse of electrolyte counterions 
onto the guest charge, the hard core of radius a, appearing in the dimen- 
sionless combination a (j3.12|) . has to be attached to the guest charge and 
only afterwards on can consider the limit a — *■ 0. Each quantity related 
to the guest particle, evaluated in the presence of the hard core, exhibits 
some relevant a-corrections with respect to its pointlike value which do not 
vanish, but on the contrary diverge, in the limit a — * 0. As Q increases, 
the hierarchy of relevant a-corrections becomes more complex and among all 
cr-corrections one is dominant in the sense that it fully determines the given 
quantity in the limit a — > 0. The main technicality of this paper consists 
in the systematic generation of the relevant a-corrections as the value of Q 
increases and in the determination of the dominant a-correction. 

The (excess) chemical potential of the guest charge, /Xgf, was the thermo- 
dynamic quantity of interest in Section 3. In the case of the pointlike guest 
charge, jJq is expressible in terms of the expectation value of the exponential 
field in the 2D sine-Gordon model, see Eq. (j2.12j) . The explicit result for 
this expectation value is available in the whole guest-charge stability region 
< (3Q < 2, see Eqs. (I2.13|) and (j2.14|) . The nontrivial extension of this 
result to the region (3Q > 2 is possible due to the existence of the reflec- 
tion relations (|3.4j) and (|3.5|) . These reflection relations exhibit an infinite 
sequence of singularities at the equidistant Q-points {Qn}^ =0 given by Eq. 
()3.7|) . The chemical potential in the presence of the hard core around the 
guest particle is the subject of Section 3.2, where the relevant a-corrections 
to the pointlike h*q are systematically generated. These a-corrections re- 
move the artificial singularities of the pointlike /i^ at the points {Q^}™ =0 , 
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which is an important consistency check of the generation procedure. The 
hierarchy of the relevant a-corrections becomes more and more complicated 
as the bare charge Q increases, however, the a-correction, dominant in the 
limit a —>■ 0, is easily detectable for arbitrary- valued Q. 

The crucial Section 4 is devoted to the evaluation of the renormalized 
charge Q ren of the guest particle. The case of the pointlike guest charge is 
briefly reviewed in Section 4.1. The inclusion of the hard core around the 
guest particle, studied in Section 4.2, requires the application of both short- 
distance (the OPE method) and large- distance (the form-factor method) 
asymptotic expansions of two-point correlation functions in the 2D sine- 
Gordon model. The systematic generation of the relevant er-corrections with 
respect to the pointlike renormalized charge is similar to the one outlined 
in Section 3.2 for the guest-charge chemical potential. The final result for 
the renormalized charge ()4.29|) is very simple. In the problematic region 
(3Q > 2, due to the condensation of electrolyte counterions onto the guest 
charge, the renormalized charge Q ren is the periodic function of the bare 
charge Q with period 1 and its value oscillates between the two extremes 
given by Eqs. (|4.3Uj) and (|4.31|) . In the Poisson-Boltzmann scaling regime, 
these two extreme values coincide and one recovers the standard Manning- 
Oosawa type of counterion condensation with the uniform saturation value 
of f3Q ren = 4/7r in the whole region j3Q > 2. From this point of view, the 
nature of the results in the PB scaling regime is fundamentally different from 
that of the exact result obtained at a strictly finite (nonzero) value of the 
inverse temperature (3. 

Appendix 

/r\\ ————— 

The integral Jq , defined in Eq. f!3.30|) , can be straightforwardly transformed 
to the form 

4 2 ) = 4vr /'drr 1 ^ 2 - 4 ^ 2 T dr' (r') 1 " 4 ^ 
J o Jo 

x dip 
Jo 

The expression in the square bracket can be expanded in the ratio (r'/r) 
by using Gegenbauer polynomials {C^(t)}, defined as the coefficients in the 
power-series expansion of the function [3*7j 

—A 00 

(l-2tx + x 2 )~ =^C J A (t)x j . (A.2) 

3=0 



r \ I r 
1 - 2 cos tp — + — 
\r \ r . 



(a.i; 



25 



Integrating then over r' and r, one ends up with 

71 ^ 1 r2n 



■I 



(2) 



l + 6 2 -2Q6 2 ^2- -j ./n 



d<pCr 2b \ CO s(p). 



(A.3) 



Gegenbauer polynomials with the cosine argument are expressible as 

n\t \ ^ r(/c + A)r(z + A) n 

(fc+i=i) 

which implies 



/•2tt 

y o d ^ c, 2j+i( cos ^) = °> 



2tt 



2tt 



dy? C^,- (cos 



The consideration of these relations in Eq. (|A.3|) leads to 



r(j + A) 
j!r(A) 



J, 



(2) 



71 



1 + 6 2 - 2Q6 2 



3=0 L 



2Qb 2 +j 



r(j - 2& 2 ) 

j!r(-26 2 ) 



(A.5) 
(A.6) 

(A.7) 



Rewriting in Eq. ()A.3j) the jth ratio 1/(2 — + j) as the integral 
Jq 1 dr r 1_4< ^ 6 the summation over j in Eq. (|A.7j) can be expressed as 

2 



2> 



J'=0 



2g& 2 + j 



r(y - 2b 2 ) 
j!r(-26 2 ) 



d 2 rr~ 4Q6 ll-r 



46 2 



7T J r <l 



(A.8) 



The integral on the rhs of Eq. ()A.8|) resembles the Dotsenko-Fateev one 
ji(—Qb 2 , b 2 , b 2 ) [see the definition ()2.2U|) ]. however, the domain of integration 
is restricted to the unit disc. It is trivial to show that 



3l (-Qb 2 ,b 2 ,b 2 " 



d 2 r 



r<l 



r -AQb 2 _j_ r 4(Qfe 2 -6 2 -l) 



1 4b 2 



(A.9) 



In the special case of equal exponents —4Qb 2 = 4(Qb 2 — b 2 — 1), i.e., when 
Q = Qc + 1/2 = Ql, one finds that 



Jl (-Q* 1 b 2 ,b 2 ,b 2 ) = 2 ( d 2 rr- 4 ^ fe2 |l-r| 4fe2 . 

Jr<l 

The comparison of this relation with Eq. (jA.8|) implies 



3=0 1 



2Q\b 2 +j 



T(j-2b 2 ) 
j\T(-2b 2 ) 



2tt 



(A.10) 



(A.ll) 
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